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基于 双向 流 固 耦合 方法 的 东风 局 特性 数值 研究 
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摘 要 风扇 转子 工作 在 高 压 比 、 高 转速 、 高 负荷 的 级 环境 中 ,实际 流动 中 存在 风扇 的 长 叶片 与 流体 间 强 烈 的 耦合 作用 . 
针对 风扇 实际 工作 中 存在 的 流 固 耦 合 问题 ， 本 文通 过 基于 时 域 推进 的 双向 流 固 磷 合 方法 对 风扇 转子 的 流 固 耘 合 问题 进行 
求解 , 并 在 此 基础 上 对 该 转子 的 超声 速 非 失速 颜 振 现象 进行 研究 . 分 析 结果 表明 , 考虑 灯 合 效应 的 转子 特性 线 整体 向 大 流 
量 方向 偏 移 0.7%， 效 率 降低 约 0.76%， 总 压 比 基本 保持 不 变 ， 稳定 工作 裕 度 略 有 减 小 ; 转子 在 设计 转速 下 近 效 率 最 大 点 
处 , 可 能 会 发 生 基于 叶片 一 阶 固有 模 态 振动 的 超声 速 非 失速 额 振 现象 . 
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Numerical Investigations on Characteristic of Fan Based on Two-Way 


Fluid-Structure Interaction Approach 


ZHANG Shuail2 GAO Li-Min'? ZHENG Tian-Longl? DENG Wei-Min!? 
(1. School of Power and Energy, Northwestern Polytechnical University, Xi'an 710072, China; 
2. Collaborative Innovation Center of Advanced Aero-Engine, Beijing 100191, China) 
Abstract There is strong coupling between the long fan blades and the fluid in actual flow field, 
because the fan rotor works in the high pressure ratio, high rotation speed and high load level 
environment. In this paper, the fluid-structure interaction(FSI) problem of the fan rotor is solved 
by two-way fluid-structure interaction method based on the time marching, and the supersonic 
unstalled flutter phenomenon of the fan rotor is studied. The numerical results show that the fan 
rotor characteristic line shift 0.796 towards the large flow direction, the efficiency is reduced 0.76%, 
the total pressure ratio remains the same and the stable operation margin decreases slightly when the 
coupling influence is considered. The supersonic unstalled flutter based on the first-order natural mode 
vibration of the blade may occur near the maximum efficiency point at the design speed of the fan rotor. 
Key words fan; fluid-structure interaction; geometric deformation; aerodynamic performance; un- 
stalled flutter 


0 引 言 


随 着 现代 先进 航空 动力 装置 性 能 的 迅猛 提升 ， 。” 题 的 研究 大 多 是 通过 流 固 耦 合 的 分 析 方 法 来 实现 的 ， 
航空 发 动机 向 着 高 推 重 比 、 高 负荷 、 高 性 能 的 方向 。” 即 采用 耦合 求解 非 定 常 气动 力 和 结构 体 瞬 时 变形 的 
不 断 发 展 。 对 于 承担 工 质 加 功 功能 的 风扇 /压气 机 ， ”方法 ,这 种 方法 能 够 提高 数值 模拟 的 精度 ， 特 别 是 
级 压 比 越 来 越 高 ， 叶 片 长 度 越 来 越 长 ， 厚 度 越 来 越 。 在 非 线性 影响 强度 较 强 的 跨 声速 区 域 ， 因 为 该 方法 


jl 


薄 ， 叶 尖 切 线 速 度 越 来 越 高 ， 非 定常 气动 负荷 不 断 能够 反映 流体 流动 与 结构 振动 的 耦合 物理 过 程 ， 不 
增 大 .叶片 在 气动 力 和 离心 力 共 同 作用 下 的 流动 不 ”但 考虑 到 结构 变形 引起 气动 力 变化 的 非 线性 ， 还 可 
稳定 问题 、 气 动弹 性 问题 等 更 加 突出 ;而 叶片 变形 ” 以 处 理 结构 动力 学 方面 的 结构 体 大 变形 、 几 何 及 材 
后 使 得 叶 型 偏离 最 佳 理 想 设 计 而 影响 风扇 /压气 机 。” 料 非 线 性 等 问题 O. 而 随 着 计算 机 技术 的 迅猛 提升 ， 


的 气动 性 能 ， 严 重 时 会 对 叶片 的 结构 强度 及 振动 计算 流体 力学 /计算 固体 力学 (CFD/CSD) 耦合 计算 
性 产生 较 大 影响 了 方法 也 得 到 了 迅速 的 发 展 B, Wilson 等 9 用 非 线 

目前 国 内 外 针对 叶轮 机 械 中 的 叶片 气动 弹性 问 性 气动 弹性 模型 研究 了 跨 声速 风扇 叶片 在 不 同 工 作 
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状态 下 ， 几 何 变形 对 性 能 的 影响 及 叶 尖 相对 安装 角 ”求解 方法 、 油 流 模型 等 的 正确 及 合理 性 ， 对 该 风月 
的 变化 , 发现 叶 片 在 不 同 转速 下 的 变形 对 气动 性 能 ”转子 的 整 机 在 设计 转速 下 进行 定常 流 场 特 性 数值 模 
影响 显著 . W. 图 2 所 示 为 该 风扇 两 级 叶片 的 流 道 及 计算 网 格 . 
随 着 计算 力学 的 发 展 ， 数 值 计 算 方法 在 叶轮 机 ”风扇 的 单 通道 四 排 叶 片 网 格 总 数 150 万 ,不 考虑 各 
械 叶片 设计 过 程 中 扮演 着 重要 的 角色 C 可 以 有 效 。 级 转 、 静 子 的 叶 顶 间隙 ， 其 进出 口 条 件 保 持 和 实验 
预测 叶片 在 实际 工作 过 程 中 存在 的 颜 振 等 气动 弹性 ”环境 条 件 基本 一 致 ,两 级 转子 的 转速 均 为 设计 转速 ， 
问题 S, 本 文采 用 基于 时 域 推进 的 双向 流 固 耦 合 数 ”进口 总 压 为 101325 Pa, 总 温 为 288.15 K, 同时 对 第 
值 计算 方法 对 NASA CR-120859 风扇 第 一 级 风扇 转 ”二 级 静 子 的 背 压 进行 调节 ; 气体 动力 学 方程 的 求解 
子 在 设计 转速 下 的 气动 弹性 问题 进行 数值 计算 研究 ， 采用 求解 定常 的 雷诺 时 均 Naiver-Stokes(N-S) 方程 
探究 风扇 转子 叶片 的 静 气 动弹 性 变形 规律 及 其 对 风 的 方法 ,清流 模型 采用 考虑 气体 前 切 应 力 的 SST i 
扁 转 子 性 能 、 内 流 场 的 影响 ; 并 在 此 基础 上 , 进一步 。 流 模型 。 
研究 该 风 遍 转 子 的 超声 速 非 失速 颤 振 现 象 ， 为 叶片 
设计 和 风扇 实际 工作 中 避免 此 类 气动 弹性 不 稳定 现 
象 提供 数值 计算 与 理论 借鉴 . 


1 计算 模型 建立 


1.1 算 例 验证 
本 文 所 采用 的 数值 模拟 理论 及 方法 均 基 于 


NASA CR-120859 风扇 ， 此 风扇 是 包括 适度 的 高 

叶 尖 速度 、 高 气动 载荷 叶片 、 多 国 弧 避 型 等 先进 设 (00 H2 两 级 风 庆 计算 风格 

计 概 念 的 两 级 风 B. 该 风 转子 的 第 一 级 转子 叶片 Fig. 2 The calculation grid of two stage fan 

的 弦 长 较 宽 、 展 弦 比 较 大 、 叶 片 较 薄 并 且 带 防 振 吓 

EE 展 纱 比 相对 PHA el 
较 小 , 该 风扇 的 各 级 动 、 更 叶片 的 实物 如 图 1 ERR. 对 比如 图 3、 图 4 所 示 . 从 对 比 图 中 发 现 , 实验 测量 


与 数值 计算 的 绝热 效率 与 总 压 比 在 变化 趋势 上 完全 
吻合 ， 在 数值 上 略 有 差异 ， 实 验 测量 的 绝热 效率 略 
高 于 数值 计算 值 ， 而 数值 计算 的 总 压 比值 略 高 于 实 
验 测量 值 . 实验 测量 的 出 口 测 点 位 于 第 二 级 静 子 出 
口 截面 ， 而 数值 计算 的 效率 计算 截面 位 于 风扇 整 机 
的 延伸 段 出 口 ， 这 导致 实验 测 基 的 静 温 值 大 于 数值 
计算 ，, 实验 测量 绝热 效率 略 大 于 数值 计算 值 , 在 几 
可 建 模 过 程 中 对 该 风扇 进行 简化 ， 数 值 计算 过 程 中 
不 考虑 风扇 转子 的 第 一 级 转子 的 叶片 防 振 突 肩 (如 


E e 
T 
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(a) First-stage blade (b) First-stage vane 图 1(a) 所 示 ) 及 各 级 转 、 静 子 的 叶 顶 间隙 等 因素 , 这 
些 均 可 能 使 该 风扇 转子 的 数值 计算 总 压 比 大 于 实验 
WEA. 


综 上 所 述 ， 本 文 所 选取 的 数值 计算 方法 、 满 流 
模型 、 转 子 几何 模型 简化 等 经 过 检验 , 风扇 转子 的 气 
动 性 能 和 实验 测量 的 变化 趋势 及 数值 基本 吻合 ， 能 
够 满足 本 文 对 于 转子 叶片 看 合 特性 研究 的 要 求 ， 
于 受 困 于 数值 计算 的 计算 量 、 转 静 交 界面 数据 传递 、 
动 欧 干涉 等 因素 ， 本 文 在 接 下 来 的 叶片 契合 特性 研 


(c) Second-stage blade (d) Second-stage vane 

M T 人 Hbk HA o B 
图 1 NASA CR-120859 RRUA ETAIT LYE 究 中 均 选 用 考虑 叶 顶 间隙 、 不 带 减 振 突 肩 的 NASA 
Fig. 1 Blades physical pictures of NASA CR-120859 CR-120859 风扇 的 第 一 级 转子 . 


1.2 本 文 计算 模型 
为 验证 本 文 所 选取 的 几何 模型 简化 方法 、 流 场 NASA CR-120859 风扇 第 一 级 转子 叶片 长 达 
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8 期 


0.236 


dl. 


TJ 


实际 转动 过 程 中 


响 变形 大 、 


m， 进 口 轮 载 比 达 到 


比 风 


P, KA 
耦合 特性 明显 ， 


张 帅 等 : 基于 双向 流 固 


.4， 叶 片 的 设计 参数 见 


ARRETA RZEK, 在 转子 


Ei 叶 


的 基础 上 而 


特性 的 影响 


究 静 气动 变 
， 选 择 该 两 级 风 


究 对 象 , 该 转子 的 三 维 几何 模型 如 图 5. 


片 受 离心 力 及 气动 力 影 
因此 为 在 流 固 耦合 分 析 
形 对 风扇 性 能 、 流 场 及 振动 
扇 的 第 一 级 转子 作为 研 
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Fig. 3 Comparison of the fan adiabatic efficiency 
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4 风扇 整 机 总 压 比 对 比 图 


EET f 


1.01 


Fig. 4 Comparison of the fan total pressure ratio 


zi 


格 生成 是 数值 计算 中 很 重要 的 环节 ， 其 质量 


决定 着 数值 计算 能 否 完成 及 其 结果 是 否 可 靠 ， 而 双 


TET RE D 


耦合 数值 计算 对 于 流体 域 及 固体 域 的 网 格 质 


合 方法 的 某 风扇 特性 数值 研究 


推进 的 非 定 常 计算 ， 计 算 域 的 网 格 是 
变化 的 动 网 格 ， 对 于 本 文 而 
转子 叶片 ， 其 叶 尖 变 
杂 ， 在 数值 计算 过 程 中 很 容易 出 现 负 
值 计算 无 法 继续 进行 。 经 过 大 量 的 网 
数量 级 网 格 的 数值 计算 发 现 , 流体 域 
下 比 较 敏 感 ， 在 数值 计算 过 程 
而 非 结构 化 网 格 的 抗 畸变 能 力 相 


变 
格 ， 


究 的 带 叶 
E RB, 


随 着 时 间 不 断 
硕 间 院 的 风 
叶 顶 间隙 网 格 复 


TR 3 


格 ， 导 致 数 
型 及 不 同 
网 格 对 于 叶片 
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Pa 


FP 很 容 
对 


易 产 生 负 网 
EER, XIF 


风扇 在 近 失 速 或 是 近 堵 塞 的 工 况 都 能 很 好 的 适应 ， 
流体 域 的 非 结 构 化 网 格 节点 数 为 621354， 单 元 数 为 


2322831， 固 壁 处 附 下 


层 层 数 为 15, 区 


体 壁 


层 网 格 尺度 满足 求解 方程 对 于 v^ 值 的 要 求 。 
表 1 NASA CR-120859 第 一 级 风 扁 转子 主要 参数 


Table 1 The main parameters of the first stage 


面 的 第 一 


fan rotor 

设计 参数 数值 

叶片 数 28 
叶 高 /m 0.2361 
设计 转速 /(r/min) 10720 
叶 尖 切线 速度 /(m/s) 441.96 

叶 顶 间隙 /mm 1.5 
设计 流量 /(kg/s) 83.55 
设计 效率 0.894 
设计 压 比 1.787 


图 


本 文 基 于 


寺 域 推进 的 双向 流 固 
解 非 定常 雷诺 平均 N-S 方程 
体 ， 瞬 态 求解 格式 是 基于 动 
后 差分 格式 , 详细 求解 器 参数 设置 如 表 2, 


5 第 一 级 转子 三 维 几 何 模型 


Fig. 5 The three dimensional geometric model of the first 


stage rotor 


IR 


上 


固体 


的 风 


基 要 求 更 高 .双向 流 固 耦 合 的 计算 方法 是 基于 时 域 


解 基于 瞬 态 动力 学 方程 


域 的 结构 瞬 态 分 析 采 


形 对 流体 域 流 场 的 影响 ， 


性 主要 关注 的 是 四 


耦合 解法 采 
， 流 体 介 质 选 取 理 想 气 
网 格 方法 的 二 阶 欧 拉 向 


男 体 


体 


因此 选取 先 求解 流体 


IR 


BK 


的 直接 积分 方法 . 在 双 
的 耦合 时 域 推进 方法 中 ， 流 体 域 的 非 定常 计算 与 


j 相 同 的 时 间 步 长 ， 而 流 
体 域 与 固体 域 的 瞬 态 求解 是 顺序 进行 的 . 本 文 研 
扇 转 子叶 片 的 耦合 特 ; 


究 
BE 


成 再 


chinaXiv:201711.02568v1 


1686 工程 热 4 
求解 固体 域 的 顺序 
表 2 求解 器 参数 设置 


Table 2 The parameter settings of solver 


计算 参数 数值 
济 流 模型 SST 
最 小 时 间 步 (t/s) 1,0x10— 
总 时 间 (T/s) 0.12 
虚拟 时 间 步 (Step) 10 
参考 压力 (Pief/Pa) 101325 
进口 总 温 (Tr1/K) 288.15 
进口 总 压 (Pri/Pa) 0 
出 口 背 压 (Ps /Pa) Variables 
2 风扇 转子 耦合 特性 分 析 
对 NASA CR-120859 风 扁 转子 进行 非 定常 计算 
和 基于 时 域 推进 的 双向 流 固 耦合 计算 ， 对 计算 得 到 


的 转子 性 能 参数 进行 时 均 化 处 理 , 图 6、 图 7 所 示 为 
两 者 的 总 性 能 参数 对 比 图 . 叶片 的 气动 弹性 变形 导 
致 转子 的 等 转速 线 发 生 明显 的 偏 移 ， 耦 合 后 的 特性 
1 线 整体 向 大 流量 方向 偏 移 0.7% 左 右 ， 变 化 趋势 旨 
本 吻合 . 效率 数值 有 所 降低 , 最 大 效率 值 降 低 0.76%, 
耦合 前 后 的 总 压 比 在 数值 上 基本 相同 ， 转 子 的 稳定 
I 


E 裕 度 略 有 减 小 . 在 双向 流 固 耦合 计算 过 程 中 ， 
于 随 着 时 间 的 推进 , 叶片 发 生 弯 扭 振动 , 导致 转子 叶 


片 的 各 截面 安装 角 变 小 , 进口 气流 角 变 大 , 流 道 的 通 
流 面积 发 生变 化 ， 使 得 风扇 转子 的 特性 曲线 发 生 整 
体 偏 移 ; 同时 , 由 于 转子 叶片 的 往复 振动 、 流 场 的 不 
断 变化 及 叶 项 间 隙 的 变化 等 因素 使 得 风扇 转子 的 等 
炉 效 率 相 比 于 耦合 之 前 有 所 降低 . 

0.94 
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Fig. 6 The entropy efficiency of fan rotor between the 


unsteady and coupling way 
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到 7 耦合 前 后 风扇 转 子 总 压 比 


Fig. 7 The total pressure ratio of fan rotor between the 


unsteady and coupling way 


2.1 流 场 分 析 

为 研究 风扇 转子 耦合 前 后 效率 及 压 比 的 变化 ， 
对 叶片 表面 的 带 压 等 参数 进行 分 析 . 图 8 给 出 
风扇 转子 叶片 在 设计 工 况 下 ， 耦 合计 算 在 叶片 
叶 高 截面 位 置 的 表面 更 压 变化 情况 ， 其 中 黑色 实 
表示 双向 流 固 耦合 的 非 定常 计算 周期 内 静 压 波动 的 
时 均值 ， 而 虚线 和 点 划 线 表示 对 应 非 定常 
片 表面 带 压 波动 的 最 大 值 与 最 小 值 . 首先 对 比分 析 
转子 叶 表 静 压 的 时 均值 ， 在 沿 着 叶 高 方向 ， 
吸 、 压 力 面 的 静 压 包 络 线 面 积 逐 渐 增 大 ， 即 从 叶 根 
到 叶 尖 ， 转 子叶 片 的 负荷 逐渐 加 大 ， 
增强 ;， 叶 表 静 压 波 动 主要 集中 于 叶片 的 吸力 面 及 靠 


近 叶 顶部 分 ， 静 压 波动 沿 着 展 向 逐渐 增强 的 同时 ， 
从 叶片 的 前 缘 波及 尾 缘 ， 而 叶片 压力 面 仅 在 叶片 前 
缘 及 98% 叶 高 的 尾 缘 处 有 较 大 的 波动 . PET ERREUR 
压 波动 变化 及 分 布 趋势 ， 与 耦合 过 程 中 叶片 的 几何 
变形 规律 紧密 相关 ， 而 叶片 变形 主要 集中 于 叶片 的 
前 缘 及 叶片 上 半 部 分 ， 


针对 转子 较 高 截面 叶 表 更 压 波动 剧烈 的 情况 ， 
如 图 9 为 时 均 化 的 叶 表 静 压 变化 趋势 对 比 图 。 图 
9(a) 中 转子 叶片 考虑 耦合 效应 的 叶 背 静 压 突变 位 置 
( 即 激 波 位 置 ) 明显 沿 弦 向 提前 7 多 左右 ， 而 压力 突 
变 的 趋势 明显 增强 ， 考 虑 耦合 效应 的 叶 背 激 波 有 所 
增强 ， 气 流 经 过 激 波 后 的 静 压 值 基本 保持 不 变 ; 图 
9(b) 中 98% 叶 高 处 的 激 波 位 置 在 叶 背 沿 弦 向 提前 了 
10% 左 右 ， 激 波 强 度 也 有 所 增加 ， 但 激 波 后 的 更 压 
值 基本 保持 不 变 ;， 而 耦合 后 的 叶片 吸 、 压 力 面 的 静 


201711.02568v1 


chinaXiv 


面积 有 所 减 小 , 转子 叶片 做 功 全 


效率 


有 


风 


AH 


以 上 分 析 ， 考虑 双向 流 固 耦合 效 
的 激 波 强度 有 所 增加 ， 激 波 位 置 沿 纺 
前 移 ， 导 致 叶片 通道 内 流动 损失 有 所 增加 ， 局 
EH IUIS, 


所 降低 ， 


3k 帅 等 : 


基于 双向 流 固 


EHA I, 


IM TESESR 


不 


时 
转 


使 得 考虑 看 合 效应 的 风月 


固体 域 计算 结果 及 分 析 
户 转 子叶 片 在 工作 过 程 中 承受 


和 气动 力 载荷, 为 了 使 得 转子 叶片 在 耦合 过 程 中 ， 


未 收敛 到 


能 小 ， 


以 及 耦合 
历程 尽 可 能 短 ， 通 过 改变 臣 


区 敛 过 程 中 叶片 几何 变 
计算 过 程 达到 稳定 收敛 位 置 的 时 让 


子 


较 大 的 离心 力 


E Et o REUS 


体 叶片 的 初始 阻尼 比 


实现 上 述 要 求 . 图 10 所 示 为 风扇 转子 叶片 在 不 同 
中 分 析 发 现 , 不 


始 阻 


尼 比 下 结构 的 响应 历程 . 从 图 


ZH 


来 
初 
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图 8 风扇 转子 在 典型 叶 高 处 表面 静 压 变化 包 线 图 


blade surface 


的 初始 阻尼 比 下 ,叶片 弹性 变形 收敛 型 


初始 阻尼 比 越 小 ， 叶 片 弹 性 变 
比 越 小 ， 
于 叶片 流体 域 网 格 质 量 要 求 越 高 ， 
越 容易 出 现 负 网 格 . 当 叶 片 的 初始 
叶片 的 静 弹 性 变形 响应 的 收敛 


形 收敛 到 稳定 幅 值 的 
时 间 历 程 越 短 ; 耦合 计算 的 收 剑 过程 中 ， 初 始 阻尼 
收敛 过 程 中 的 不 稳定 振动 


u 


敛 过 程 中 的 网 格 最 大 变形 量 远 高 于 


The static pressure change package of the typical 


相同 的 幅 值 ; 


幅 值 就 越 大 ， 大 
在 计算 过 程 中 就 
阻尼 比 为 0.1 时 ， 
最 短 , 但 其 在 收 
其 他 ， 这 对 叶片 


流体 域 的 网 格 质 量 提出 了 很 高 的 要 求 。 因此， 需要 


选择 一 个 兼顾 耦合 过 程 
变形 量 的 叶片 初始 阻 
图 11 为 初始 阻尼 


比 e=0.3 


4 时 间 
2 IA Mi (e= 


叶 尖 截面 前 缘 、 中 部 及 尾 缘 (LE, Mid 


历程 与 网 格 最 大 
3). 
T, X 


记 转 子叶 
及 TE) 
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在 离心 力 载 荷 与 气动 力 载荷 作用 下 的 响应 。 从 变 ”叶片 前 缘 变 形 幅 值 来 看 ， 叶 片 在 近 失 速 的 工 况 点 变 
形 幅 值 来 分 析 ， 叶 人 尖 前 缘 的 变形 幅 值 最 大 ， 叶 尖 ERRAK. 在 转子 的 等 转速 工作 线 上 , 叶片 叶 尖 截面 
中 部 变形 幅 值 次 之 ， 叶 尖 尾 缘 的 变形 量 最 小 ， 可 最 大 变形 量 随 着 流量 的 增加 而 逐渐 减 小 ,特别 是 叶 
以 发 现 叶 片 的 振动 形式 是 扭转 ， 从 振动 的 稳定 性 来 。” 人 尖 前 缘 变 形 量变 化 更 加 剧烈 ,这 说 明 此 处 离心 力 载 
分 析 ， 叶 尖 截 面 尾 缘 最 先进 入 稳定 状况 ， 而 变形 ” 和 荷 与 气动 力 载 荷 在 叶 尖 作用 方向 相同 ， 气 动力 引起 
量 最 大 的 叶 尖 截面 前 缘 振 动 幅 值 较 大 并 且 收 敛 速 度 。 的 叶片 变形 量 随 流量 的 增加 而 逐渐 减 小 .而 在 风月 
gH. 转子 的 近 堵 塞 工 况 时 ， 离 心力 引起 的 叶片 变形 量 占 
总 变 的 比重 急剧 上 升 ， 甚 至 在 某 些 工 况 出 现 叶 

T€ - 片 的 ， F 仅 离心 力 载 荷 作 用 的 总 变形 量 ， 
uL. Unsteady | 这 说 明 此 时 有 可 能 出 现 离心 力 与 气动 力 在 叶 尖 截面 
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Fig. 9 The relative pressure at the typical blade surface 


distributes along the direction of the normalized chord length 
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作用 相反 的 状况 ， 气 动力 的 作用 减弱 了 离心 力 对 于 Da(Inter-Blade Phase Angle, IBPA) 变化 图 .叶片 
叶片 的 作用 。 在 一 阶 模 态 振动 时 ， 气 动 阻尼 系数 在 IBPA-334.3* 
附近 出 现 负 值 ， 即 不 稳定 区 域 ， 而 第 二 、 三 阶 模 态 
0.008 T 对 应 的 气动 阻尼 系数 均 为 正 值 . 通过 计算 分 析 发 现 ， 
E —L— Mi 对 于 该 风扇 转 子叶 片 , 一 阶 模 态 为 其 最 危险 的 模 态 ， 
E zd 且 其 最 危险 的 IBPA 大 概 在 334.3° 附近 . 通过 控制 
Eoo IBPA=334.3° 不 变 , 改变 风扇 转 于 叶片 的 振动 频率 ， 
M 得 到 其 不 同 叶片 振动 频率 对 应 的 气动 阻尼 系数 ， 通 
$ 上 过 插值 得 到 转子 气动 阻尼 系数 为 零 的 非 失 速 颜 振 临 
Soon. 界 点 。 
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图 12 叶 尖 监测 点 总 变形 量 随 风扇 转子 工 况 变化 


Fig. 12 The total deformation of the blade tip monitoring 


point changes with the operating condition of fan 
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由 文献 [9] 知 叶轮 机 械 的 主要 颤 振 形式 有 : K 
EMR, AAMIR, GARRE, A100 型 
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58—9& Ug ETI ETER RRENA. 由 
文献 [9],[ 
子 效率 较 高 的 工 况 ， 选 取 该 风扇 转子 
的 近 效 率 最 高 点 处 的 工 况 进行 研究 . 轮机 械 
子叶 片 以 往 的 实验 研究 发 现 ， 转 子叶 片 颤 振 受 其 前 
三 阶 振动 模 态 影响 较 大 ， 而 高 阶 振动 模 态 由 于 本 身 
对 应 的 叶片 振动 频率 很 高 ， 而 实际 中 叶片 发 


In 
HT 


的 频率 较 小 且 振动 幅 值 较 小 中, 所 以 文章 主要 基于 
能 量 法 D 对 该 风扇 转子 的 前 三 阶 模 态 的 振动 现象 
进行 着 重 分 析 . 

通过 对 风扇 第 一 级 转子 的 前 三 阶 模 态 的 非 定常 


气动 力 做 功 及 气动 阻尼 系数 计算 来 达到 分 析 转 子叶 
片 稳定 性 问题 . 模 态 分 析 得 到 该 转子 叶片 固有 的 一 
阶 模 态 频率 为 189.34 Hz, 二 阶 模 态 频 率 为 599.32 Hz, 


三 阶 模 态 频率 为 967.69 Hz. 图 13~15 分 别 为 风 局 
转子 前 三 阶 模 态 对 应 的 气动 阻尼 系数 随 叶 间 相 位 
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到 13 f—189.34 Hz 一 阶 模 态 下 叶片 气动 阻尼 


Fig. 13 Aerodynamic damping of blade in first-order mode 
when f—189.34 Hz 
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图 14 f=599.32 Hz 二 阶 模 态 下 叶片 气动 阻尼 
Fig. 14 Aerodynamic damping of blade in second-order mode 
when f=599.32 Hz 
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的 形式 发 生 超声 速 非 失 速 颤 振 。 而 通过 改变 固体 
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到 15 f=967.69 Hz 三 阶 模 态 下 叶片 气动 阻尼 
Fig. 15 Aerodynamic damping of blade in third-order mode 
when f—967.69 Hz 
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图 16 IBPA—334.3* 一 阶 模 态 气动 阻尼 系数 


Fig. 16 Aerodynamic damping of blade in first-order mode 
when IBPA=334.3° 


依据 线性 插值 对 气动 阻尼 系数 零点 进行 插值 ， 
当 气 动 阻尼 系数 为 零 时 ， 其 对 应 的 叶片 振动 频率 
为 f=201.34 Hz( 如 图 16 所 示 )， 即 为 一 阶 模 态 对 
应 的 颤 振 边 界 。 对 于 该 风扇 第 一 级 转子 而 言 ， 第 
一 阶 模 态 是 其 前 三 阶 模 态 中 最 不 稳定 的 模 态 ， 所 
以 可 以 认为 一 阶 模 态 额 振 边界 即 为 该 风扇 转子 的 非 
REMA. 以 上 分 析 表 明 ， 转 子 在 设计 转速 
下 近 效 率 最 大 点 处 ， 一 阶 模 态 频率 低 于 f=201.34 
Hz 时 ， 转 子叶 片 可 能 会 以 一 阶 模 态 、IBPA=334.3? 


叶片 


的 阻尼 (增设 减 振 突 肩 等 ) 可 以 使 叶片 避免 


非 失速 颤 振 的 发 生 ， 以 达到 使 风扇 稳定 工作 的 
Hif. 


425 论 


通过 基于 时 域 推进 的 双向 流 固 耦 合 方法 对 


NASA CR-120859 第 一 级 风扇 转子 进行 非 定常 数 
值 模拟 ， 并 在 此 基础 上 对 该 转子 的 超声 速 非 失 速 
颤 振 特 性 进行 了 研究 ,通过 计算 分 析 得 到 以 下 主要 


N 


阶 模 态 频率 低 于 f=201.34 Hz 时 ， 转 子叶 片 可 能 
会 以 一 阶 模 态 、IBPA=334.3" 的 形式 发 生 超声 速 
JEREMI. 通过 改变 叶片 的 结构 阻尼 (增设 减 
振 突 肩 等 ) 或 材料 阻尼 可 以 避免 其 非 失 速 颜 振 的 
AE. 


向 大 流量 方向 整体 偏 移 OT CAU E He uy dS 
效率 有 所 下 降 ， 最 大 效率 点 降低 0.76%; 转子 的 
总 压 比 基本 保持 不 变 ; 转子 的 稳定 工作 裕 度 有 所 


片 吸 、 压 力 面 的 静 压 包 络 线 围 成 面积 有 所 减 小 ， 叶 
片 做 功能 力 减弱 ， 


1) 考虑 双向 流 固 耦 合 效 应 的 风扇 转子 性 能 曲 


2) 考虑 耦合 效应 的 计算 结果 表明 ， 叶 背 处 激 波 


位 置 向 叶片 前 缘 移动 ， 最 大 位 移 量 达到 1096, mj 
波 的 强度 也 略 有 增强 ， 叶 片 通道 损失 增加 ; 转子 叶 


3) 转子 在 设计 转速 下 的 近 效 率 最 大 点 处 ， 一 
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